#include <stdio.h>
#include <omp.h>
#include <math.h>
#define  N 1000000000


int main(int argc, char *argv[])
{
double a = 0;
double b = 1.0;
double h = (1.0/N);
int i=0;
double z = 0;
double result = 0;
// #pragma omp parallel for private(z) shared(h) reduction(+:result)
for(i=0;i<N;i++)
{
	z = i*h;
result += (z*sin(z*z))*h;


}

printf("result: %f \n",result);

return 0;
}
